High-density K nuclear systems with isovector deformation 
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Using a phenomenological KN potential which reproduces A(1405) as an / = bound state of 
KN, we investigated deeply bound kaonic nuclei, ppnK~ and 8 BeK~, with the method of Antisym- 
metrized Molecular Dynamics. Our calculations show that strongly bound kaonic nuclear systems 
with unusual exotic structures are formed around the K~, which attracts the surrounding nucleons 
to an extremely high-density assembly and induces a proton-neutron separation, "isovector defor- 
mation" . 
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Recently, exotic nuclear systems involving a K (K~ 
or K°) as a constituent have been studied theoretically 
by Akaishi and Yamazaki ]jj (to be abbreviated as AY), 
who constructed KN interactions phenomenologically so 
as to reproduce low energy KN scattering data (2| , kaonic 
hydrogen atom data || and the binding energy and de- 
cay width of A(1405), which is asserted to be an I = 
quasi-bound state of KN. These interactions are charac- 
terized by a strongly attractive 1 = part. Such proper- 
ties of the AY KN interactions are consistent with those 
obtained by using a chiral SU(3) effective Lagrangian 
by Weise el ai. ||. With these interactions AY calcu- 
lated the bound systems of ppnK~, ppnnK~ and 8 BeK~, 
showing the following characteristics: i) The 1 = KN 
interaction is strong enough to shrink the nucleus against 
the nuclear imcompressibility. ii) The binding energies 
are extremely large due to the strong 1 = KN in- 
teraction helped by the nuclear shrinkage effect so that 
the bound states lie below the threshold of the main de- 
cay channel to £7r, thus inferring the presence of quasi- 
stable discrete bound states (width (Tk) < binding en- 
ergy (-Bk))- Furthermore, Yamazaki and Akaishi H have 
shown that the (K~,7r _ ) reaction as a source of A* can 
produce various exotic K bound states in non-existing 
proton-rich nuclei, such as p 2 ( 2 He), p 3 ( 3 Li), p 3 n ( 4 Li), 
p 4 n 2 ( 6 Be) and p 5 n 4 ( 9 B). Since the presence of a K - 
in nuclei is expected to cause a drastic change in the 
nuclear systems, it is vitally important to perform the- 
oretical studies of such exotic systems without any con- 
straint from the existing common knowledge of nuclear 
physics (shell and cluster structure, nuclear shape, nu- 
clear density, proton-neutron distributions, etc.). Thus, 
we started a series of calculations, employing the method 
of Antisymmetrized Molecular Dynamics (AMD) || 
based on the AY KN interaction. The present paper re- 
ports the first results on ppnK~ and 8 BeK~. 

The AMD method has a great advantage: it treats "a 
nucleus + K~" as a full (^4 + l)-body system without 



any assumption concerning a priori structure (shell or 
clucter) . Actually, AMD has succeeded to explain a vari- 
ety of exotic shapes and clusters in light unstable nuclei 
0. If a K nuclear system prefers a structure where the 
proton distribution differs from the neutron distribution, 
such a structure will be dynamically formed in the AMD 
treatment. 

In the usual AMD, a single nucleon wave function, \(pi), 
is represented by a single Gaussian wave packet. In the 
present study, however, it is given as a superposition of 
several Gaussian wave packets as follows, in order to de- 
scribe the system more precisely, and the K is represented 
in the same way: 
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where \(JiTi) is a spin-isospin wave function, which is |p ] 
),|p !),|n |), or |n |); |tk) indicates the isospin wave 
function of K. We use two Gaussian wave packets for a 
single nucleon (tin = 2) and five for a K (n^ = 5). We 
must antisymmetrize the nucleon wave functions, which 
are then combined with a K wave function as 



det[\<pi(j))] <g> \<pk) 
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Finally, we project the intrinsic wave function |$) onto 
the eigen-state of parity: |$ ± ) = |$) ±V\<&). We utilize 
|"I> ± ) as a trial function, which contains complex vari- 
ational parameters {C^,Z l a , C^,Z^} and a real varia- 
tional parameter, v. They are determined by the fric- 
tional cooling equation with constraint 
The Hamiltonian, 



H = T + y NN + Vk N + Vc - T G , 



(4) 
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(c) Proton (d) Neutron (e) K" 

FIG. 1: Calculated density contours of ppnK~. Comparison between (a) usual 3 He and (b) 3 HeK~ is shown in the size of 5 by 
5 fm. Individual contributions of (c) proton, (d) neutron and (e) K~ are given in the size of 3 by 3 fm. 



is composed of the kinetic energy, T, two-body NN in- 
teraction, Vnn, KN interaction, Vkn, and Coulomb force, 
Vc- The center-of-mass motion energy, Tq, is subtracted. 
In the kinetic energy and the cm. motion energy, we 
treat the mass difference between a nuclcon and K _ ex- 
actly. Since the system is likely to be extremely dense 
due to the strongly attractive KN interaction, we can use 
no existing effective interactions. Therefore, we use the 
Tamagaki potential (OPEG) || as a bare NN interaction, 
and the AY KN interaction as a bare KN interaction, and 
construct the effective NN-central force and the KN force 
with the G-matrix method Since the Tamagaki po- 
tential reproduces the NN phase shifts up to 660 MeV in 
laboratory energy we expect that it is applicable to 
extremely high-density states of relevant concern. 

We have calculated systems of ppnK - and 8 BeK~. 
These systems were previously treated with a Brueckner- 
Hartree-Fock method by AY ]jj . We summarize our cal- 
culations for ppnK~ and 8 BeK~ in Table I. We also show 
the AMD results of the normal 3 He and 8 Be for a com- 
parison. 3 He was calculated with the Volkov No.l force 
H and 8 Be with the MV1 Case-3 force 0. However, 
we did not perform an angular momentum projection for 
simplicity. Tk is the width for decaying to A-7T and E-7T, 
and is evaluated by calculating the expectation value of 
the imaginary potential contained in the effective AY KN 
interaction with the wave function obtained by the AMD 
calculation. "ppnK - t" an( f " 8 BeK~ t" are the results 
of AY. 

Our result shows that the ppnK" (T = 0) system is 



TABLE I: Summary of the present calculations. p(0): nu- 
cleon density at the center of the system. Runs'- root-mean- 
square radius of the nucleon system, v: width parameter of 
a Gaussian wave packet used in the calculation. j3: deforma- 
tion parameter for the nucleon system. ppnK - ' and 8 BeK~ 
. AY's results. 
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3 He 


5.95 




0.14 


1.59 


0.20 


0.0 


ppnK - 


113 


24 


1.39 


0.72 


1.12 


0.19 


ppnK - t 


116 


20 


1.10 


0.97 






8 Be 


48.72 




0.10 


2.46 


0.205 


0.63 


8 BeK~ 


159 


43 


0.76 


1.42 


0.52 


0.55 


8 BeK- t 


168 


38 


~ 0.85 









bound by 113 MeV (B K =105 MeV below the 3 He+K~ 
threshold). Since the obtained state lies under the T,tt 
threshold (—80 MeV), the T,tt channel, which is the main 
decay-channel, is closed and the decay width, Ppc, be- 
comes narrow to be 24 MeV. The present result is very 
similar to the AY prediction: Bk = 108 MeV and 
Tk = 20 MeV. We have not considered the decay width 
from the non-mesonic decay (KNN — > AN/EN), but ac- 
cording to AY it is estimated to be about 12 MeV pj. 
The width of ppnK - remains still narrower than that of 
A(1405), even when the non-mesonic decay is taken into 
account. 
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(c) Proton (d) Neutron (e) K" 

FIG. 2: Calculated density contours of 8 BeK~. Comparison of the density distributions of (a) usual 8 Be and (b) 8 BeK _ is 
shown in the size of 7 by 7 fm. Individual contributions of (c) proton, (d) neutron and (e) K~ are given in the size of 4 by 4 
fm. 



Surprisingly, the central density ( "uncorrelated den- 
sity") of the system amounts to 8.2-times the normal 
density due to the shrinkage effect. Figs. 1(a) and (b) 
show a comparison between 3 He and 3 HeK _ . In order 
to see how the bound K changes the nucleus in more de- 
tail we show the calculated density distributions of the 
constituents in Figs.l(c)-(e). Apparently, the proton dis- 
tribution is more compact than the neutron distribution. 
This phenomenon is attributed to the property of the KN 
interaction. Table II shows how protons and a neutron in 
ppnK~ contribute to the kinetic energy and the expecta- 
tion value of the KN interaction, and also to each root- 
mean-square radius. This table together with Fig.l can 
be interpreted as follows. Since the K~p interaction is 
much stronger than the K _ n one, the protons distribute 
compactly near the K _ . Though the kinetic energy of 
the protons increases, the total energy decreases due to 
the strongly attractive K _ p interaction. On the other 
hand, the neutron is not favoured by shrinking, since the 
K~n interaction is not attractive enough to compensate 
for the increase of the kinetic energy. 

Now we proceed to 8 BcK~. As expected, 8 Be is ex- 
tremely shrunk due to the presence of a K~. The K~ is 
bound by 159 MeV (B K = 104 MeV with respect to 8 Be). 
The binding energy of K - and the decay width obtained 
by our calculation are very similar to those by AY, who 
employed the aaK~ model fjj]. We compare the den- 
sity distribution of 8 BeK~ with that of the normal 8 Be 



TABLE II: Contributions of proton and neutron in ppnK - . 
(T): kineic energy per nucleon. {Vkn}'- expectation value 
of the KN interaction per nucleon. -R rms : root-mean-square 
radii of the proton and neutron distributions. 





(T) 


(Vkn) 


-Rrms 




[MeV/u] 


[MeV/u] 


[fm] 


proton 


69.4 


-170.2 


0.70 


neutron 


59.0 


-42.4 


0.75 



in Figs. 2(a) and (b). 8 BeK~ becomes much more com- 
pact than 8 Be, and attains a high density, 4.5-times the 
normal nuclear density. As shown in Table I, the root- 
mean-square radius decreases dramatically from 2.46 fm 
to 1.42 fm. Nevertheless, an a-cluster like structure still 
remains in 8 BeK~. We can see that the relative distance 
between the a clusters becomes short and that the a clus- 
ter, itself, in 8 Be is also shrunk. The shrinkage of each a 
cluster is indicated by that of the Gaussian wave packet, 
as shown in Table I. The width parameter, v, is changed 
from 0.205 fm -2 for the normal 8 Be to 0.52 fm~ 2 for 
8 BcK~. The deformation parameter, (3, becomes slightly 
smaller. 

Figs.2(c)-(e) display contour mapping of the calculated 
density distributions of 8 BeK~ in more detail. Although 
the two-a like structure, as shown in (b), persists, it is 
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opments of nuclear physics will follow. This new field of 
"bound-K nuclear spectroscopy" will create a new dense 
and cold environment, in contrast to hot and dense 
states expected to be realized by heavy-ion collisions. 
Restoration of chiral symmetry in a dense nuclear 
medium might be observed for the KN system in the 
same way as in deeply bound pionic states jl2|, O, fT!] . 
The bound-K nuclear spectroscopy will yield unique 
information about a transition from the hadron phase to 
a hitherto untouched "quark phase" . It will also provide 
experimental grounds for possible kaon condensation 
and strange matter formation p5|. 



FIG. 3: Calculated density distributions of proton, neutron 
and K~ in 8 BeK~ along the horizontal axis in Fig. 2. For 
comparison the density distribution of s Be is also displayed 
as " 8 Be". 

found that the proton distribution (c) is significantly dif- 
ferent from the neutron distribution (d). As in the case of 
ppnK~, the protons feel a strongly attractive force from 
the K~, which stays in the center of the total system 
(Fig. 2(e)), and are attracted nearer to the K~ than the 
neutrons. Namely the proton distribution is separated 
from the neutron distribution in spite of N = Z. A more 
quantitative presentation of the proton and neutron den- 
sity distributions is given in Fig. 3. The K produces a 
new type of neutron halo; it is peripherally distributed, 
yet with quite a high density. We name such a new type 
of deformation as an "isovector deformation" . It should 
be noted that AMD enables us to predict such an isovec- 
tor deformation without any assumption about resultant 
structure. Such a deformation cannot be seen in the cal- 
culation by AY because they assumed the existence of a 
clusters. 

In summary, we have shown from our AMD cal- 
culations that a K~ forms a deeply bound nuclear 
state by producing an extremely high-density nuclear 
system with isovector deformation. No such nuclear 
systems at "zero temperature and high density" have 
been created so far. Whereas an experiment to identify 
3 HeK~ from 4 He(stopped K~, n) jl], is under way 
at KEK, a new possibility to make use of (K~,7r~) 
reactions to populate proton-rich systems has been 
proposed ||. We are extending our AMD calculations 
to such more exotic systems, as AMD is the only 
possible method to investigate their dynamics without 
any a priori structural assumption. If a deeply bound 
kaonic nucleus is confirmed experimentally, new devel- 
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